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INTRODUCTION 


Three or more isotopic elements as strontium and lead are widely utilized in the dating 
of rocks. Thus a need arises for faster data processing with statistical methods. The usual 
way of meeting this demand in the case of new instruments is an incorporation of a com- 
puter, via an interface anda direct statistical adjustment of the data obtained during mass 
spectrometric measurements. However, many older instruments still working in scientific 
laboratories are not equipped with such devices and statistical processing of the data 
obtained must be performed in a computer centre. The purpose of this paper is to eliminate 
such inconvenient step by using a simple programmable calculator which can be even 
communicated with the mass spectrometer employing a home-made interface similar to 
that described by Hatas and Skorzynski [2]. 2 

The Authors have accomplished this goal using inexpensive TI Programmable 59 calcu- 
lator produced by Texas Instruments. This calculator has up to 100 registers availble for 
data storage or up to 960 for program storage. The memory area can be optionally divided 
between these two needs. The calculator reads and writes magnetic cards. 


TERMINOLOGY AND METHODS 


Let us consider a mass spectrum of an element for which three peaks are measured. As 
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an example may serve the mass spectrum of strontium (Fig. 1) in which the isotopic ratio 
8751/85 5г is of geological interest. 

However. if the measurements will be limited merely to this ratio, the obtained result 
can be alter-ed by both systematic and random errors. In order to reduce the systema- 
tic error the final value of isotopic ratio is normalized by the following equation [3]: 

US: Е SeS (ES/AS 8 ' 
85 Sr (855/050), : Р (0 


where index mm refers to measured ratios and п to a normalized one (know with a higher 
accuracy). In order to reduce the random errors, the measured ratios, being utilized in 
the above equation, should be adjusted by the Leas Squares Method. For this reason all 
major peaks should be measured and three ratios must be calculated. 

To unify the notation for stronAium and lead isotopes (Fig. 1, 2), the heights of peaks 
86 Sr ®7Sr and 88 $г will be denoted asx;, хз and x4, and their ratios as: 


Ja = X3/x2 
Js = ха/х. . (2) 
Je = Xa/Xs 


These three values are not independent, but are connected by the following rigorous 
equation: 


Nns — ns = 0, (3) 


where n's denote the true values of respective y's. 

Each of the determined ratios will be considered as a sum of true value and of an 
attached error. A set of the measured values, generally, does not satisfy the rigorous 
equation. According to the Method of Least Squares, however, the weighted sum of 
squares of errors can be minimalized, thus providing the most reliable estimates of the 
detrmined values which satisfy the rigorous equations (see, for instance [1]). 

In the formulation of general equations it is convenient to use the ordinary matrix 
notation. In this notation the determined values, in our example the strontium isotope 
ratios, will be denoted as one column matrix, i.e. vector: 


У4 ВЯ 
у=| уз = | 851/557 |. (4) 
Ye 88 Sr/87 Sr 


The variances of y 5, i.e. squared standard errors, will be represented by diagonal ele- 
ments of the variance-covariance matrix: 
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Fig. 1. Mass spectrum of strontium; the heighst of peaks are proportional to 

the abundances of the respective isotopes. Three major peaks are measured 

to reduce statical error of %51/% $ ratio; note that the haight of the **Sr 
peak has been reduced by one-tenth 


Pp 


0 || 


Fig. 2. Mass spectrum of lead, note that the height of the 
2% Pb has been enlarged by tenfold 
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о 


оо 
Ct =| 0 о oj- (5) 
0 оо 
The first estimates of vector % (representing true values) are calculated from the follo- 
wing general formula derived by the use of Lagrange’s multiplers method ([1], chapter 9): 


y = у —C,BG'c (6) 


where c is the ‘vector of discrepancy’, the elements of which are obtained by replacing 
7s in the rigorous equations by the respective values of y's, 


С = BC, B, (7) 


and B is the matrix the rows of which are derivatives of each of the rigorous equations 
with respect to ту; values taken at ту = yj. Prim in the above formulas indicates that а 
given matrix is transposed. 

As usual in nonlinear problems the first estimates may not sufficiently satisfy the rigo- 
rous equation, therefore iterative procedure should be used: the second approximation of 
the vector g can be calculated again from the above equations by putting y instead of the 
initial vector y in matrix B, vector c and equation (6). Further approximations, if 
necessary, can be calculated in the same way. 

Finally, the variance- covariance matrix of improved ў vector can be calculated by 
applying to equation (6) the propagation of errors formula ([1], chapter 4—5): 


Cy = Cy — CyBG" BC, . (8) 
Square roots of diagonal elements in this matrix are standard errors of improved 
values y , whereas symmetric nondiagonal elements are convariances between the individual 


pairs of the improved values. From the element of matrix Су one may calculate the esti- 
mates of correlation coefficients: 


ew Iu ams al (9) 


In the case of strontium isotope ratios the above defined matrices and formulas of the 
Method of Least Squares may be rewritten as follows: 

B = Qs. FL У): (10) 

ВС, = (0455, —03, 08), (11) 


С = 04$ + 0$ + 08 » (12) 
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УУ У peine. 
~ 6” Js 
y — y = —а? ? (13) 
G 2 
0674 | : 
4972.-999, OORT aye | 
C, = Cy = с —0401Y, , oj , —00у, ы (14) 


04037455 , —O8OR Ys, 081 


Here matrix С is reduced to a number, therefore С! equals simply 1/G. Note that 
matrix in equation (14) is symmetric, hence only three non-diagonal elements must be 
calculated. 

Example 1. Let us adjust a set of measured strontium ratios represented by the following 
vector у and variance-co-variance matrix: 


Уг | ' 07122] 0.81X 10$ 0 0 2 
у = |у; | = | 8.3625 |, Cy = 0 4.84X 10 5 0 * 
Ув. 11.7316. 0 0 7.84Х 1076 


The initial ‘discrepancy vector’ is —0.00725. Using equations (12), (13) and (14) the 
first and the second approximations were computed. After the first approximation the 
‘vector of discrepancy’ was 1.9X 107? while after the second one it became as small as 


—2.7X 107", so the calculations were stopped and variance-covariance matrix of improved 
values was calculated: 


"Totis 0.06 0.38 —0.44 
y-|8362|, Су = 10% | 038 4.65 022|- 
11.7319. | —0.44 0.22 7.58. 


Actually, in this case, one step of approximation is enough since the differences between 
the same elements of vector y, after the first and the second approximations, are smaller 
than 10°’. Note that the adjustment process essentially reduces the variance of geochemi- 
cally important ratio уд = ® Sr/%* Sr. 


ADJUSTMENT OF LEAD ISOTOPE RATIOS 


In the case of lead spectrum an example of which is shown in Fig. 2, the accurate deter- 
mination of three isotope ratios, namely ?Pp/?9* Pb, 297 pp/?9* Pb апа 205 pp/29* Pb, is 
important in many geo- and cosmochemical studies. These ratios are well known as X, Y 
and Z values in scientific papers. In our notation they will be represented by уу, Y2 and 
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y3 values, respectively. Further three ratios which can be determined from four peaks at 
the mass spectrum, are usually chosen in the following order: 


Ja Sei, р — СЫРЫ, ie = 0200 p 20 prp 


We will denote the heights of peaks by x4, x2, хз and x4. Let vector y represent six 
ratios to be determined: 


У Х2/х! 
Age x3/x, 
у = || = |xahy] * (15) 
Уз хз/х 
Js Xa/X4 
Je Xa/xs | 


From the above matrix equation it follows that these measurements are connected by 
three independent equations. Further calculations will be based on the following set of 
rigorous equations: 


1114 — 7 = 0 
nins — пз = 0 (16) 
1216 — Тз = 0 


This set provides the design matrix 


Jae SNe PRS Nye. 0 0 
0 Ve ale 0 0 Уг 
and the `уесюг of discrepancy’ 
iX. — У2 
CEN S33! ve (18) 
УЕ УЗ 


which can be used in the data adjustment with general formulas (6) and (8). Let us consi- 
der a numerical example. 

Example 2. For a lead sample the following data were experimentally obtained, 
represented by the vector 


16.07 
15.49 
у = 136.78 
0.9651 
2.2891 
2.3726. 
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and variance-covariance matrix 


Cy = 1074 


The initial ‘discrepancy vector’ is 


0.0192 
c= 0.0058 | * 


.—0.0284. 


while after the first step of iterative process (equation (6)) its elements are as small as 
about 10-7. After the second step they are comparable with machine error = 107 . Cal- 
culations are, therefore, stopped at the second step of approximation, and the improved 
values and their variance-covariance matrix are obtained: 


16.0624 211 200 476 —0.23  —0.49 0.06 
15.5002 200 197 458 0.24  -0.02 -—0.62 
Уна 36-2005 |6; Су = 107*| 476 458 1109 —0.08 123 1.47 
0.9650 —0.23 0.24 -0.08 0.029 0.028 —0.042 |* 
2.2893 —0.49 —0.02 1.23 0.028 0.146 0.082 
2 3108 0.06 —0.62 1.47 0.042 0.082 0.189 


This example illustrates that estimates of covariances between the improved 1 ratios ys 
Уз, Ye and the others are very small. In contrast to these values, the ratios J, , Y2, Уз are 
strongly correlated with correlation coefficients approaching the limit value p = 1. 
Such a behaviour of the elements of variance-covariance matrix and the'strong convergence 
of iterative process indicate that our problem is nearly linear. Actually, as the height 
of the minor peak is by one order smaller than the rest, its small variation produces large 
and strongly correlated variations in Jı š Va and Ўз ratios. Moreover, these three ratios 
show percentage errors by one order higher than three others. Hence 4, ns and mg in the set 
of rogorous equations (16) may be considered as nearly constant. 

One may conclude from this discussion that the problem of the adjustment of lead 
isotope ratios can be divided into two smaller problems: adjustment of y4, ys, yg and, then, 
У, Уз, Уз ratios. The first problem is identical with that demonstrated above for strontium 
isotopes, whereas the second problem may be solved in the same way. However, it includes 
two independent rigorous equations and three measurements to be adjusted. Hence, the 
method of elements ([1], chapter 9—3) is more convenient than the method of Lagrange 
multiplers used above. 

Putting the adjusted values Y4, Уз, Ув into equations (16), one may obtain three 
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linear equations, but only two linearly independent. Let us choose the following set for 
further calculations: 


туа -m = 0 


= (19) 
uyd?c FA 
Taking 71 as the element, from equation (19) we obtain 
71. rd " 
=] = |у | т = Ет. (20) 
У; 


ЦЕ] 


To solve our adjusting problem one column matrix Е, thus obtained, must be put in the 
following expressions: 


m 1 ~ ~ 
FF’ i is Ja в 4 
Су = ni^ SH 92. * "3. (21) 
s Jays У а 
where 
1 y yn 
Н = ЕСЕ =—+— +25 , (22) 
1 2 oj 
and 
A y/o} 
у = C6 у = Су | 2/03 |. (23) 
¥3/03 
Note that all correlation coefficients followed from the expression‘s (21) equal 


unity, as expected. 

Example 3. In order to compare the results obtained in this very simplified way with 
those obtained by the method demonstrated above, let us consider again the same set of 
initial values y and Cy as in Example 2. Our simplified approach provides the following 
first estimates: 


16.0624 209 202 478 0 0 0 
15.5002 202 195 461 0 0 0 
У = | 36.7715 |, Су = 10* |478 461 1095 0 0 0 
0.9650 0 0 0 0.029 0.009 —0.043 
2.2893 0 0 0 0.029 0.147 0.083 


2137291 no 0 0  —0.043 0.083 0.190_ 
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This example clearly shows that in the case of lead isotopes, and also in other cases when 
one peak is small in comparison with the other tobe measured, both methods provide nearly 
identical adjusted ratios. Also the standard errors and correlation coefficients, derived 
from simplified Су, are close to those derived from a strict variance-covariance matrix. 
Even those correlation coefficients which correspond to covariances are estimated here 
as zero. Therefore, a calculator program will be presented for the simplified procedure 
only. The program for the strict method on TI Programmable 59 calculator is available 
from the authors. 


THE CALCULATOR PROGRAM 


The program listed in the Appendix has been written according to Owner’s Manual 
for TI Programmable 58/59. It consists of 316 steps and uses 33 registers, therefore it can 
be written on one magnetic card without the reset of the normal partitioning of memory 
storage area in TI Programmable 59. The flow diagram of the program is shown in Fig. 3. The 
program consists of two basic parts. The first part (steps 000 through 127) beginning with 
label A enables the adjusting of *"Sr/9*Sr, 85Sr/86Sr, Sr/®7Sr or 207pp/206 pp. 208pp206pp. 
208 pp/297 Pb ratios, і.е. уг, ys, yg in our unified notation. 

The second part (steps 128 through 240) provides the first estimates of 206 Pb/™ Pp, 
207 pp/?9* Ph, 208 pp/?95 Pb ratios noted here as y, , уз, ys values. To this part belong three 
subroutines (labels B, B' and E) including instructions for printing improved values with 
their standard errors. Subroutine C is designed for a convenient i ntroduction of initial 
data: ratios and their standard errors. Tables 1 and 2 cotain user instructions and register 
contents. Figure 4 shows two examples of standard printouts. 

As numerous free steps are available at the normal partitioning of storage area, we have 
extended our program to a subroutine (not listed here) for calculation of initial ratios and 
errors directly from the heights of peaks. This capability can be utilized in an instantaneous 
adjustment of readouts of digital voltmeter, being proportional to heights of peaks, if 
the calculator is interfaced with a mass spectrometer. Such automatization has been 
accomplished in our laboratory for the measurements of isotope ratios in the case of light 
elements and we are going to extend this technique for geochemically important heavy 
elements. 
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MAIN PROGRAM 


LblA 
Calculate matrix BCy, CyB'and G“ 
based on equations (10), (11) and (12) 


Multiply matrix CyB'by G^' Master 
Library Рет 03E 


Calculate Y,, ¥,, У, using equation (13) 


057 


058 


086 


087 


Calculate Cy using eq. (14) Pgm 03C 
used 


126 


127 
Calculate matrix FF and Н”! based on 
eqs. (20), (21), (22) 174 
175 
Calculate Cy (equation (21)) Рет 03E 
weet 186 
187 
Calculate y (equation (23)) Рет 03С” 
= 216 
217 
Printing instructions, subroutines B, В” 
and E used 239 


Fig. 3. The flow diagram of the program for adjusting strontium and lead isotope ratios 


SUBROUTINES 


240 


Lbi C 
Data introducti 
ata introduction 275 


LbIE 
Round-off instructions 


276 


283 


284 


| ГЫ В | 
Print vect 
rint vector y 298 


ГЫ В” 
Calculate errors and print 


299 


316 
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Tab. 1. User Instructions 


Step Procedure Enter Press Display 
1 Initialize (© 28 
2 Enter elements of vector у У R/S У, 


y, R/S Ya 
y,  R/S ys 
MR R;S Ya 
У; R/S Ys 


У R/S 6 
3 Enter standard errors 0, R/S о? 
0; R/S о? 


vc, R/S o) 
в. R/S о} 
0; R/S о? 
2, R/S o; 
4 Perform the adjustment of the entered values A  seenotel 


NOTES: 1. Subroutine A provides standard printout, shown in Fig. 4, with 4 decimal places. One 
may reset this number occupving step 279 of the program. 

2. For strontium isotopes perform all above steps with any number as y, , Y3, y, and в,, 05,0; 
values. 

3. Elements of variance-covariance matrix may be recalled after computations from appropriate 
registers, see Tab. 2. 


Tab. 2. Register Contents 


Ro Ra 

Ка Ry | 

Raz R,, } Used 
Ros 7 Ri; 

В | Used в. 

Ке | R,, Used 
Ros | Ris 

R,, Pointer R4 

Ry Used , Cn К, У, 

R, Used ,c4 R4 У, 

R, Used , Cu Ry Ys 

R, m7 s Ča Ray». 
R,, ar re Cn К, Ws, У: 
Ris Ce Сва: En R4, Yas Уз 
Ry, » Св Ry c, 

R, Used ‚еш К, с, 

Ry Ас _ 
R,, Used „у, Rs; си, Cas 
Ris Used » У: К, Css, Css 
К, Cs „Ӱз Ras Cone Сев 


NOTES: 1. A fraction of memory registers are used twice: the first symbols represent number 
stored when the first part of the program is running while the second symbols represent numbers 
which can be recalled after run of the second part. 

2. Symbols су were used for the elements of variance-covariance matrix of improved values. 
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28. 0000 


16. 0700 У, 
ae ed 15. 4900 У, 
Sr rem 36. 7800 y, 
onip: 0. 9651 у 
.3625 Ys 2.2891 VQ 
11.7316 У 2.3726 y* 
. c < 6 
"7 ED 0. 2498 0, 
36. 0. 2358 ©; 
1296. 0000 x 0. 5892 © 
A ao =) 0. 0021 e? 
. (eins no M 
: S 0.0045 Б 
п. 0028 ci 0. 0055 6, 
25. 1604 16.0624 Ў, 
20. 0435 15. 5002 $, 
235. 1457 м 36.7715 
D. 7128 ; 0.9650 $ 
8. 3622 2. 2893 
9 > : 
11.7319 $ 2. 3723 9, 
r X 
5 9537 D. 1445 & 
4. 243 0. 1395.6 
49. 7561 м n, 3303 = 
0. 0002 2 0.0017 # 
о. 0022 Dg 0. воза O 
. 0028 5, 0. 0044 6, 


Еф. 4. Standard format of the printed results; the time of calculation is about 1 min 
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049 
050 
051 
052 
053 
054 
055 
056 
057 
058 
059 
060 
061 
062 
063 
064 
065 
066 
067 
068 
069 
070 
071 
072 
073 
074 
075 
076 
077 
078 
079 
080 
081 
082 
083 
084 
085 
086 
087 
088 
089 
090 
091 
092 
093 
094 
095 
096 
097 


01 1 

42 STO 
04 04 
00 0 

42 STO 
06 06 
36 PGM 
03 03 
15 E 

43 RCL 
31 31 
65 x 

43 RCL 
33 33 
75 - 

43 RCL 
92-32 
95 = 

42 STO 
15 15 

36 PGM 
03 03 
BBC 
43 RCL 
EI 
22 INV 
44 SUM 
31 31 
43 RCL 
12 12 
22 INV 
44 SUM 
32 32 
43 RCL 
13 13 
22 INV 
44 SUM 
33 33 
43 RCL 
16 16 
42 STO 
№5 15 
36 PGM 
03 03 
18 C 
43 RCL 
11 11 
22 INV 
44 SUM 


140 
141 
142 
143 
144 
145 
146 


APPENDIX 
PROGRAM LISTING 
37 37 14785 + 
43 RCL 148 43 RCL 
17 17 149 36 36 
42STO 150 35 ИХ 
15 15 151 65 x 
36 PGM 152 43 RCL 
03 03 15332 32 
Lome 154 42 STO 
43 RCL 155 10 10 
12 12 156 42 STO 
22INV 157 14 14 
44 SUM 158 33 X? 
38 38 159 42STO 
43 RCL 160 16 16 
11 11 16195 = 
42STO 162 35 ИХ 
19 19 163 42 STO 
43 RCL 164 05 05 
18 18 165 43 RCL 
42STO 166 31 31 
i545. че ба б5 вх 
36 РСМ 168 43 RCL 
03 03 169 32 32 
182€ 170 95 - 
43 RCL 171 42 STO 
1359131 017295515 
221NV . 173 42S STO 
44 SUM 174 13 13 
39 39 17503 3 
01 1 176 42 STO 
42STO 177 04 04 
08 08 178 04 4 
43 RCL 179 42 STO 
34 34 180 25 25 
351/X 18100 0 
85 + 182 42 STO 
43 RCL 183 06 06 
35 35 184 36 PGM 
351/X 185 03 03 
65 x 186 15 E 
43 RCL 187 43 RCL 
31 31 188 34 34 
42STO 189 35 ИХ 
09 09 19065 x 
42STO 191 43 RCL 
11 11 19228 28 
33 X? 19395 = 
42 STO 194 42 STO 
12 12 19521 21 


196 43 RCL 
197 35 35 
198 35 1/X 
199 65 x 
200 43 RCL 
201 29 29 
202 95 = 
203 42 STO 
204 22 22 
205 43 RCL 
206 36 36 
207 53 ИХ 
208 65 x 
209 43 RCL 
210 30 30 
211 95 = 
212 42 STO 
213 23 23 
214 36 РСМ 
215 03 03 
216 18 С 
217 01 1 
218 07 7 
219 42 STO 
220 06 06 
221 12 B 
222 03 3 
223 01 1 
224 42 STO 
225 06 06 
226 12 B 
227 98 ADV 
228 08 8 
229 42 STO 
230 06 06 
231 17 B. 
232 03 3 
233 07 7 
234 42 STO 
235 06 06 
236 01 1 
237 42 STO 
238 25 25 
239 17 В’ 
240 76 LBL 
241 13 C 
242 06 6 
243 42 STO 
244 07 07 


245 02 2 
246 08 8 
247 42 STO 
248 06 06 
249 91 R/S 
250 72 ST* 
251 06 06 
252 15 E 
253 69 OP 
254 26 26 
255 97 DSZ 
256 07 07 
251102 102 
258 49 49 
259 06 6 
260 42 STO 
261 07 07 
262 98 ADV 
263 91 R/S 
264 15 Е 
265 33 X? 
266 72 ST* 
267 06 06 
268 69 ОР 
269 26 26 
270 97 052 
271 07 07 
272 02 02 
273 63 63 
274 98 ADV 
275 92 RTN 
276 76 LBL 
277 15 Е 
278 58 FIX 
279 04 04 
280 99 PRT 
281 22 INV 
282 58 FIX 
283 92 RTN 
284 76 LBL 
285 12 B 
286 03 3 
28742 STO 
288 07 07 
289 73 RC* 
290 06 06 
291 15 E 
292 69 OP 
293 26 26 


88 Stanistaw Hatas, Piotr Staniaszek 


294 97 052 299 76LBL 304 73 RC* 309 25 25 314 03 03 
295 07 07 300 17 B' 305 06 06 310 44 SUM 315 04 04 
296 02 02 301 03 3 306 34 TX 311 06 06 316 92 ВТМ 
297 89 89 302 42STO 30715 Е 312 97 052 
298 92 КТМ 303 07 07 308 43 RCL 313 07 07 


STRESZCZENIE 


Praca zawiera prosty i praktyczny algorytm pozwalający poprawiać stosunki izotopowe dla pier- 
wiastków złożonych z kilku (3—5) izotopów oraz umożliwiający znaczną redukcję błędów losowych 
obliczanych w ten sposób stosunków. 

Algorytm wykorzystuje metodę najmniejszych kwadratów i ma charakter iteracyjny. W przedsta- 
wionych przykładach wystarczyły dwa kroki iteracyjne do tego, by różnice między kolejnymi przy- 
blizeniami poprawionych elementów stały się porównywalne z błędem maszynowym (ok. 1071). 
Algorytm został zrealizowany na programowanym kalkulatorze kieszonkowym TI Programmable 59. 
Opracowany program jest szczególnie użyteczny przy poprawianiu stosunków izotopowych strontu 
i ołowiu. 

Dodatki zawierają sieć działań programu, przykładowe tabulogramy wyników oraz wydruk 
programu. 


РЕЗЮМЕ 


В работе излагается простой алгоритм, который позваляет значительно поправлять 3Kcne- 
риментальные данные по изотопным отношениям элементов составленных на 3—5 изотопов. 
В алгоритме, который имеет итерационный характер, используется метод найменьших квадра- 
тов. В Передставленных вычислениях достигзется разницу между последовательными при- 
ближениями поправленых отношений порядка 10 !! (машинной ошибки) в двух итерацион- 
ных шагах. 

Алгоритм реализован на программированном калькуляторе ТИ-59. Разработанная програм- 
ма особенно полезна для поправляния изотопных отношений стронта и свинца, В приложениях 
представлены: диаграмма действий программа и полная форма программа. 


Ztozono w Redakcji 2 ХІ 1981 roku. 


